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The crossing of a transition state in a multidimensional reactive system is mediated by invariant geometric 
objects in phase space: An invariant hyper-sphere that represents the transition state itself and invariant 
hyper-cylinders that channel the system towards and away from the transition state. The existence of these 
structures can only be guaranteed if the invariant hyper-sphere is normally hyperbolic, i.e., the dynamics 
within the transition state is not too strongly chaotic. We study the dynamics within the transition state for 
the hydrogen exchange reaction in three degrees of freedom. As the energy increases, the dynamics within 
the transition state becomes increasingly chaotic. We find that the transition state first looses and then, 
surprisingly, regains its normal hyperbolicity. The important phase space structures of transition state theory 
will therefore exist at most energies above the threshold. 
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I. INTRODUCTION 

Transition state theory 1-4 (TST) is the cornerstone of 
reaction rate theory. It assumes that reactant and prod- 
uct regions in phase space can be separated by a dividing 
surface that all reactive trajectories must cross once and 
only once. If this condition is satisfied, TST allows one 
to calculate the (classical) reaction rate exactly. Other- 
wise, it provides an upper bound to the rate. For this 
reason, great effort has been devoted to the construction 
to a surface that is recrossing free or at least minimizes 
recrossings. 

Pollak and Pechukas 5-9 identified the optimal dividing 
surface in collinear exchange reactions with two degrees 
of freedom: It is given by the projection into configu- 
ration space of an unstable periodic orbit and is called 
a periodic orbit dividing surface (PODS). This surface 
will yield an exact reaction rate as long as there is only a 
single periodic orbit in the transition region 7 . Otherwise, 
the family of periodic orbits can be used to compute both 
upper and lower bounds to the reaction rate 8 . The phase 
space structures in the transition region that lead to the 
failure of TST have been studied in detail for a variety of 
two-dimensional reactive systems (see, e.g., Refs. 10-18). 

In systems with more than two degrees of freedom, 
a recrossing-free dividing surface was found only much 
later 19 . Such a surface exists only in phase space, not in 
configuration space. It is bounded by a high-dimensional 
invariant hyper-sphere that plays the role of the periodic 
orbit in the two-dimensional setting. At sufficiently low 
energies, this hyper-sphere is a normally hyperbolic in- 
variant manifold (NHIM), i.e., the dynamical expansion 
and contraction rates transverse to the hyper-sphere are 
larger then those in directions parallel to it. There are 
two important consequences of normal hyperbolicity 20,21 : 
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First, the invariant hyper-sphere will persist under per- 
turbations of the dynamical system, for example, changes 
in energy. Second, the hyper-sphere possesses stable and 
unstable invariant manifolds. These manifolds separate 
reactive from non-reactive trajectories in phase space. 
They act as reaction channels that guide the system from 
the reactant configuration towards the transition state 
and on into the product region. Their knowledge allows 
a detailed description of the reaction dynamics that goes 
well beyond a rate calculation 22-28 . 

The invariant hyper-sphere and its stable and unsta- 
ble manifolds will persist as long as the sphere is nor- 
mally hyperbolic. This can be guaranteed for energies 
just above the reaction threshold. The reaction dynam- 
ics within the transition state region is then accurately 
described by a harmonic approximation 19 . The dynamics 
within the invariant hyper-sphere is therefore completely 
regular, and the condition of normal hyperbolicity is sat- 
isfied. At higher energies, the dynamics within the sphere 
will become partially chaotic, and a breakdown of nor- 
mal hyperbolicity may result. Such a scenario has indeed 
been described in 29-31 for a model reaction. The authors 
analyse the dynamics with the help of normal form trans- 
formations. Because this procedure does in general not 
converge, it can become difficult, in particular at higher 
energies, to distinguish the properties of the underlying 
dynamical system from artefacts of the normal form. In 
this paper, we will investigate the dynamics within the 
transition state of a physical system and present a de- 
tailed description of those features that lead to a break- 
down and, surprisingly, to a subsequent reestablishmcnt 
of normal hyperbolicity. 

The invariant hyper-sphere that embodies the transi- 
tion state has customarily 19,22 ' 23,26-28 been called "the 
NHIM." Because we are interested in situations in which 
the sphere fails to be normally hyperbolic, we will avoid 
that term and call this object the central sphere. 

We study the H + H 2 exchange reaction that has also 
served as the prototypical example for the analysis of 
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transition state structures in two degrees of freedom. We 
will focus on the dynamics within the central sphere, viz., 
within the centre manifold of the saddle point that me- 
diates the exchange. The centre manifold forms an in- 
variant subsystem with two degrees of freedom. It can 
be analyzed by means of a Poincare surface of section. 

There are two fundamental periodic orbits within the 
centre manifold: a symmetric stretch periodic orbit 
(SSPO) and a bending periodic orbit (BPO). In addition, 
a secondary symmetric stretch periodic orbit (ScPO), 
that is generated by a bifurcation, plays an important 
role because it gives rise to a large regular island. These 
periodic orbits undergo a sequence of bifurcations in 
which they successively lose and regain stability. 

To decide at what energies the central sphere is nor- 
mally hyperbolic, we compute Lyapunov exponents in 
directions parallel and perpendicular to the sphere for 
both periodic and non-periodic orbits. We find that the 
SSPO is the only orbit that violates the condition of nor- 
mal hyperbolicity, and that only in a small energy inter- 
val. Nevertheless, the dynamics of the exchange becomes 
more and more complex as the energy is increased, as is 
evident already in the collinear subsystem. 

The contents of this paper are as follows. In Sec. II 
we give an overview of the low-energy phase space struc- 
tures near a saddle point. In Sec. Ill, we present the 
H + H2 exchange reaction in three degrees of freedom. 
In Sec. IV we investigate the dynamics within centre 
manifold. We describe the fundamental periodic orbits 
and their bifurcations, and we illustrate the dynamics by 
means of Poincare surface of section plots. In Sec. V we 
compute Lyapunov exponents of trajectories within the 
centre manifold and identify energy intervals in which it 
is normally hyperbolic. 



and VH = (|f , ||, If, ||). They read explic- 
itly 

Qj = P:i > Pi = - u faj for j = 1, . . . , n - 1, 

q n =Pn, Pn = A 2 g„. (5) 

The eigenvalues of the matrix associated to the lin- 
earized Hamiltonian vector field around the saddle point 



are ±A and ±io;,- where j = 1, . 



1. The pair real 



eigenvalues ±A describe the hyperbolic direction of the 
saddle point while the complex eigenvalues describe the 
elliptic directions of the saddle point, i.e., oscillations 
transverse to the reaction coordinate. 

The dynamics described by the Hamiltonian (1) has 
a stationary point at pi = qi = at energy zero. We 
will study the dynamics at a fixed energy h > above 
the reaction threshold. The energy surface is (2n — 1) 
dimensional and is given by 
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From (6), we can see the that the section through the 
energy surface at fixed q n is a (2n — 2) sphere with radius 

yh + ^-q n - Thus the energy surface is a hyper-cylinder 
S 2n - 2 x R. 

The centre manifold of the equilibrium point contains 
all trajectories that remain trapped close to the equi- 
librium for all time in the infinite future and the infinite 
past. It is given by q n = p n = 0. This surface is invariant 
because the equations of motion (5) imply q n = p n = 0. 
It has dimension 2n — 2. For a fixed energy ft > it 
intersects the energy shell in a surface that satisfies 



II. THE PHASE SPACE STRUCTURE NEAR A SADDLE 
POINT 

We now present the phase space structure of a linear 
Hamiltonian system near an equilibrium point of centre 
x centre X...X saddle type 19 . Close to the saddle point, 
the dynamics is well described by the harmonic Hamil- 
tonian 

j=i j=i 

The corresponding equations of motion are given by 



dH dH 
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dp, 1 dqj ' 

or, in terms of the phase space vector x = 
(q 1 ,...,q n ,p 1 ,...,p n ) G M 2n , by 

x=-J-VH, (3) 

where 
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This equation describes an (2n — 3) dimensional hyper- 
sphere S^ n ~ 3 . This is the central sphere (called the 
NHIM in earlier studies) that forms the bottleneck for 
phase space transport from reactants to products. It has 
stable and unstable manifolds attached to it. These are 
(2n - 2) dimensional manifolds, denoted by VF s (Sf l "~ 3 ) 
and W u (Sf l n ~ 3 ), respectively. They are given by 

71—I n—1 
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These manifolds are referred to as reaction cylinders. 
Their structure is S* 2 ™ -3 x BL The stable and unsta- 
ble manifolds have two branches, forward and backward 
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cylinders, denoted by and W b ' , respectively. 
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A recrossing-free dividing surface is a (2n — 2) dimen- 
sional hyper-sphere that is defined by setting q n — 0. It 
has codimension one in the energy shell and separates re- 
actant from product regions. Each trajectory that crosses 
the dividing surfaces crosses only once, going from reac- 
tants to products if p n > or from products to reactants 
if p n < 0. The only exceptions are trajectories within 
the central sphere, which have p n = and remain in the 
dividing surface for all times. The central sphere is an 
equator of the dividing surface and splits the dividing sur- 
face into two hemispheres with p n > and p n < that 
mediate forward and backward reactions, respectively. 

The internal dynamics of the central sphere, according 
to (5), is described in this approximation by a multidi- 
mensional harmonic oscillator and is therefore completely 
regular. In the transverse direction the central sphere is 
unstable because it is balanced near the top of an ener- 
getic barrier. It is therefore normally hyperbolic. This 
feature guarantees that both the central sphere and its 
stable and unstable manifolds persist in the full anhar- 
monic system, at least at low energies where the central 
sphere lies close to the equilibrium point and the anhar- 
monic terms are small. At higher energies, the invariant 
manifolds, assuming they exist, can be approximated via 
normal form transformations 19,23 ' 26 ' 29 ' 30 . In the present 
work, we will avoid normal forms and investigate the per- 
sistence of the invariant manifolds by direct numerical 
simulation. 



collinear configurations, which are invariant under rota- 
tions around the axis on which the atoms lie 32 ' 33 . 

The origin of this singularity can be illustrated with 

0,the help of Fig. 1. To obtain the configuration space of 
the vibrational dynamics, we have to identify all configu- 
rations of the reactive complex that can be transformed 

0,into each other by translations or rigid rotations. The 
shape of the complex can then be described by the three 
coordinates r, x and y, where r is the distance between 

0,111 and H2, y is the perpendicular distance from H3 to 
the distance r, and x is the distance from the midpoint 
of HI and H2 to the end of the line through H3 perpen- 

O.dicular to r. However, not all those configurations are 
different: (r, x, y) can be transformed (r,x,—y) by a ro- 
tation around the axis through HI and H2. To resolve 
this ambiguity, the configuration space must be restricted 
to the half space (r, x, y > 0). It has a boundary that is 
formed by the collinear configurations with y = 0. The 
dynamics must necessarily be singular at these configu- 
rations. 

We are mainly interested in studying the dynamics in 
the vicinity of the saddle point that marks the transi- 
tion region for the exchange reaction. Unfortunately, the 
activated complex is collinear at the saddle point, and 
the ensuing singularity makes it difficult to analyze the 
dynamics. To circumvent this difficulty, we regard con- 
figurations with positive and negative values of y as dif- 
ferent. This convention, which has also been employed 
in previous studies 26 ' , can physically be interpreted as 
constraining the three atoms to move in a plane. The full 
system then has six degrees of freedom, two of which cor- 
respond to the centre of mass motion and one to planar 
rotations. The remaining three degrees of freedom, which 
can be described, for example, by the three coordinates 
r, x and y of Fig. 1, describe the vibrational dynam- 
ics of the complex. Because the collinear configurations 
are not invariant under rotations in the plane, or equiva- 
lently, because configurations with positive and negative 
y cannot be transformed into each other through planar 
rotations, the symmetry-reduced phase space does not 
have singularities. It is well suited to an investigation of 
the dynamics near the saddle point. 

We will study the vibrational dynamics of the complex 
at zero angular momentum. The Hamiltonian is then 
given by 



III. THE HYDROGEN EXCHANGE REACTION 

The hydrogen exchange reaction H + H 2 — > H 2 + H 
involves three atoms. Consequently, if the atoms are 
assumed to move in three-dimensional space, the reac- 
tion is described by nine degrees of freedom. Three of 
these, which represent the centre of mass motion, can 
be separated directly. Of the remaining six degrees of 
freedom, three describe spatial rotations of the complex 
and three describe vibrations. However, the attempt to 
separate rotational from vibrational degrees of freedom 
leads to a vibrational phase space that is singular for all 
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P 2 r + j(pl+P 2 v ) + 



(xp y -yp x f 



■V(r,x,y), 
( 10 ) 

where vtlh is the mass of the hydrogen atom (tojj = 
1.00794 amu). The expression for the kinetic energy is de- 
rived from that given, for example, by Waalkens et al 2e . 
This transformation is shown in the appendix. We use 
the potential energy surface V(r,x,y) derived by Porter 
and Karplus 35 . Distances will be measured in atomic 
units (a.u.) and energies in electron volts (eV), with the 
potential energy of three isolated hydrogen atoms chosen 
as zero. 
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The potential energy V(r, x, y) has two reflection sym- 
metries: x i — y —x and y n- — y. That these transfor- 
mations must leave the potential invariant is clear from 
Fig. 1 because the three atoms are identical. The two 
reflections of configuration space are extended to phase 
space by the canonical transformations 

Px ■■ {r,x,y,p r ,p x ,p y ) i-> (r,-x,y,p r ,-p X) p y ), 

P y : (r,x,y,p r ,p x ,p y ) (r,x,-y,p r ,p x ,-p y ). (If) 

Both of these are symmetries of the Hamiltonian (10), as 
is their composition 

P x oP y : (r,x,y,p r ,p x ,p v ) h-> (r,-x,-y,p r ,-p x ,~p y ). 

Corresponding to the two reflection symmetries (11) 
there are two subsystems with two degrees of freedom. 
They contain all configurations that are invariant under 
one of the reflections. 

The reactive complex is invariant under P y if y = p y — 
0. These are precisely the collinear configurations. The 
Hamiltonian of the collinear case is given by 

H = J- [p 2 r + lp x ] +V(r,x). (12) 
ran 

Numerous researchers including Pollak and co-workers 5-8 
and most recently Inarrea et a/ 36 have studied the 
collinear hydrogen exchange reaction. 

The subsystem invariant under the reflection P x con- 
tains all axially symmetric configurations with x — p x — 
0. The dynamics within this subsystem is described by 
the Hamiltonian 

H= — [p 2 r + lpl]+V(r,y). (13) 

The saddle point of the Porter-Karplus potential en- 
ergy surface is located at the symmetric collinear config- 
uration (r,x,y) = (rs, 0, 0) with r§ = 3.40166 a.u. To 
obtain a harmonic approximation of the dynamics close 
to the saddle point, we expand the Hamiltonian (10) in 
a Taylor series up to second order. The last term in the 
kinetic energy will not contribute because it is of fourth 
order. Due to its symmetries the expansion of the po- 
tential energy cannot contain any terms of odd order in 
either x or y. Up to an additive constant the harmonic 
Hamiltonian must therefore be of the form 

H 2 = — [pl + l(p 2 x +p 2 y )]+a(r-r s ) 2 +bx 2 +cy 2 (14) 
ran 

with constants a, b, c that cannot be determined from 
symmetry considerations. Thus, the dynamics in r, x 
and y will decouple in the harmonic approximation. Be- 
cause the expansion point is a saddle, the dynamics must 
be unstable in one of the three coordinates, namely the 
reaction coordinate. In the reactant and product states 
the middle atom H3 is bound to either HI or H2, whereas 
the third atom is far away. It is therefore plausible to 
identify the reaction coordinate with the coordinate x 




FIG. 1. The coordinates for H + H2 exchange reaction. 

that brings H3 closer to one or the other atom. Indeed, 
the expansion of the Porter-Karplus potential shows that 
the coefficient b is negative whereas a and c are positive. 

As a consequence, the symmetric subsystem x = p x — 
in which the motion in the reaction coordinate is sup- 
pressed forms the centre manifold of the transition state, 
i.e. it contains all configurations in which the system os- 
cillates around the unstable equilibrium point. The sym- 
metry of the system makes it easy to identify the centre 
manifold without laborious calculations. It allows us to 
avoid the normal form calculations that are required in 
reactive systems without this symmetry. 19 ' 29,30 



IV. DYNAMICS WITHIN THE CENTRE MANIFOLD 

The central sphere that controls transport through the 
transition state at low energies can be identified with 
the energy shell within the centre manifold, as described 
in Sec. II. As we aim to investigate the breakdown of 
the low-energy phase space structures, we will start by 
studying the dynamics within the centre manifold. 

Fig. 2 shows a contour plot of the potential energy for 
symmetric configurations (i.e. within the centre mani- 
fold, with x = 0). The saddle point (sp) of the three- 
dimensional system appears as a minimum. It lies at 
(3.40166 a.u., 0) in r, y coordinates and has energy value 
—4.3504 eV. The second prominent feature of the poten- 
tial is a conical intersection ridge. It occurs at equilateral 
configurations, where y = ^-r. For these configurations 
the too lowest electronic states are degenerate. As a con- 
sequence, the potential energy surface, which gives the 
energy of the lowest state, is not smooth at the intersec- 
tion. The lowest point on the ridge occurs at r = 1.90352 
a.u., y = ±1.64849 a.u. with energy E cA , = -1.9514 eV. 
Above this energy, a new reaction channel opens in which 
the central atom (iJ3 in Fig. 1) can escape across the 
ridge, leaving the two outer atoms bound as a molecule. 
The transition across the conical intersection cannot be 
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FIG. 2. The contour plot of the potential energy surface of 
the centre manifold (x = 0) of the equilibrium saddle point 
of the Hamiltonian flow, (sp) refers to saddle point and (c.i.) 
refers to the conical intersection. 



described by classical mechanics. We will restrict our 
following investigations to energies below E c ,i_. As we 
will see, complicated dynamics develop well below this 
threshold. 

For energies between the saddle point and the con- 
ical intersection, the contour line of the potential en- 
ergy is topologically a circle. The energy shell in phase 
space (within the centre manifold) has therefore the same 
topology as it has in the harmonic approximation, i.e., it 
is a three-dimensional hyper-sphere that we have called 
the central sphere. At low energies, it is normally hy- 
perbolic. As the energy increases, normal hyperbolicity 
might, and indeed will, be destroyed. We know from 
these simple considerations, however, that the central 
sphere will persist even at energies where it is not nor- 
mally hyperbolic. 

At energies close to the saddle point where the har- 
monic approximation is accurate, the dynamics within 
the central manifold can be described by two normal 
mode vibrations, a symmetric stretch and a bend of the 
activated complex. Their frequencies can be obtained 
from a second-order Taylor series expansion of the po- 
tential, i.e., from the constants a and c in Eq. (14), as 
usspo = 4.1121 x 10 14 s" 1 and lo B po = 1-8458 x 10 14 
s _1 . Both normal mode periodic orbits are stable with 
respect to a perturbation of initial conditions within the 
centre manifold. As the energy increases, they undergo 
a sequence of bifurcations in which they lose their sta- 
bility and give rise to further stable periodic orbits, as 
illustrated schematically in Fig. 3. 

Fig. 3 shows the bifurcation diagram of the two fun- 
damental periodic orbits, the symmetric stretch (SSPO) 
and the bend (BPO) within the centre manifold. The 
first bifurcation occurs in the SSPO atBw -4.32547 eV, 
just above the saddle point energy. The SSPO undergoes 



a period doubling bifurcation: It becomes unstable and a 
new stable periodic orbit with twice the period appears. 
We will see that this periodic orbit plays an important 
role in structuring the dynamics within the centre man- 
ifold. We will call it the secondary symmetric stretch 
periodic orbit (ScPO). 

The configuration space projections of the fundamental 
periodic orbits are shown in Fig. 4 for energies E = —4.3 
eV and E — —4.0 eV. It can be clearly seen that even 
though the ScPO is generated by a bifurcation from 
the SSPO, it takes on pronounced bending character at 
higher energies. All three periodic orbits are invariant 
under the reflection P y . The SSPO is located within the 
collincar subsystem, which means that each point on the 
SSPO is invariant under reflection. This is not true for 
the ScPO and BPO. These periodic orbits are invariant in 
the sense that any point on one of these orbits is mapped 
under reflection to a different point on the same orbit. 
Periodic orbits of this type can undergo symmetry break- 
ing bifurcations that do not exist in systems without re- 
flection symmetries 37 ' 38 : A stable periodic orbit that is 
invariant under reflection turns unstable and gives rise to 
two stable periodic orbits that are not invariant, but are 
mirror images of each other. The asymmetric periodic 
orbits have roughly the same period as the symmetric 
one. 

A symmetry breaking bifurcation of the ScPO occurs 
at the energy E rj —2.5 eV. Fig. 5 shows the configu- 
ration space projections of the ScPO and the two new 
periodic orbits for energy E = —2.3 eV. It is obvious 
from the figure that the satellite orbits have lost their re- 
flection symmetry. At a higher energy E —2.2 eV the 
asymmetric periodic orbits collapse onto the ScPO again 
and the ScPO regains stability in an inverse symmetry 
breaking bifurcation. 

In a similar scenario, the BPO undergoes a symmetry 
breaking birucation at E ss —3.8 eV, and the two asym- 
metric periodic orbits thus generated collapse onto the 
BPO again and disappear at E w —2.4 eV in an inverse 
symmetry breaking bifurcation. These three periodic or- 
bits are shown in Fig. 6. The two asymmetric orbits have 
the same projection into configuration space, but, as the 
phase space figures show, they are traversed in different 
directions. 

In order to investigate the dynamics within the centre 
manifold in more detail, we choose a suitable Poincare 
surface of section. Since the centre manifold is four di- 
mensional, the surface of section will have two dimensions 
and will be easy to visualize. We pick the surface of sec- 
tion y = and use the canonically conjugate variables 
r and p r as coordinates in the surface. The remaining 
momentum p y can be determined from the energy con- 
servation condition 

H(p r ,p y ,r,y = 0) = E. (15) 

We will always choose 

Py = Py{Pr, T, E) > 
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FIG. 3. Bifurcation diagram for the system within centre manifold. The solid and dashed black lines denote stable and unstable 
periodic orbits, respectively. PD refers to periodic doubling bifurcation and SB refers to symmetry breaking bifurcation. 



at the initial point, and in computing the Poincare map 
we only consider intersections of a trajectory with the 
surface of section that havep y > 0. The SSPO lies within 
the surface of section. Indeed, it bounds the area that 
is energetically accessible at a given energy. In contrast, 
the BPO appears as the central point in the low energy 
surface of section. For various energies, the surface of 
section is shown in Fig. 7. At low energies, the inter- 
sections of a single trajectory with the surface of section 
lie on a closed curve, indicating quasi-periodic motion in 
accordance with the Kolmogorov-Arnold-Moser (KAM) 
theorem 39 . 

Some of the bifurcations of the fundamental periodic 
orbits, such as the loss and return of stability of the BPO, 
can also be seen in Fig. 7. The bifurcation of the SSPO 
has an unusual appearance because the SSPO forms the 
boundary of the surface of section. As a consequence, 
the ScPO appears at the boundary and moves toward the 
centre of the surface of section. The Poincare plots show 
only a single periodic point corresponding to the ScPO, 
as the chosen surface of section is p y > 0. A second peri- 
odic point is located in the surface p y < 0. Both periodic 
points can be seen in Fig. 8, which shows the Poincare 
surface of section r = rg = 3.40166 a. u. The SSPO in- 
tersects this surface transversely, and its bifurcations are 
therefore shown more clearly. The neighborhood of the 
SSPO in Fig. 8 has the appearance one would expect close 



to a period doubling bifurcation. Note, however, that the 
situation is different from that shown in Fig. 7(e) in the 
neighbourhood of the BPO. As indicated by the colors, 
the two periodic points appearing there belong to two dif- 
ferent periodic orbits, each of which has approximately 
the same period as the BPO. The two periodic points 
close to the SSPO in Fig. 8 lie on a single periodic orbit 
of twice the period. 

In addition to the fundamental periodic orbits, the sur- 
face of section plots show many other, longer periodic 
orbits that are not included in Fig. 3. Of these there 
are, of course, infinitely many. Most important for our 
purposes is the observation that regions of chaotic dy- 
namics appear and grow as the energy is increased. If 
the dynamics within the central sphere is chaotic, the 
central sphere might fail to be normally hyperbolic. We 
will investigate this question in the following section. 



V. BREAKDOWN OF NORMAL HYPERBOLICITY 

For the hydrogen exchange reaction, we have seen that 
we can guarantee the existence of the central sphere for 
energies up to the conical intersection ridge without hav- 
ing to rely on its normal hypcrbolicity. The full geometric 
structure of TST, however, also requires the existence of 
the reaction tubes, i.e., the stable and unstable mani- 



2.5 3.0 3.5 4.0 4.5 5.0 

r(a.u.) 



(b) 




2.5 3.0 3.5 4.0 4.5 5.0 



r(a.u.) 

FIG. 4. Periodic orbits on the potential energy of centre man- 
ifold at energies (a) E = -4.3 eV, (b) E = -4.0 eV. 

folds of the central sphere. This can only be guaranteed 
if the central sphere is normally hyperbolic. We will now 
investigate the energy range in which this is the case. 

We first need to state the condition of normal hy- 
perbolicity more precisely. It is expressed in terms of 
Lyapunov exponents, which measure the rates at which 
nearby trajectories diverge under the dynamics: The 
Lyapunov exponent associated with variations of the ini- 
tial conditions within the centre manifold (\i n t) should 
be less than the Lyapunov exponent away form the cen- 
tre manifold (X ex t)- Thus normal hyperbolicity survives 
as long as \ ex t > Xint- For small enough energies above 
the saddle energy, the internal Lyapunov exponents are 
zero because the dynamics within the central sphere is 
completely regular. 

To compute a Lyapunov exponent for an arbitrary tra- 
jectory, consider a trajectory x(t) and a neighboring tra- 
jectory x(t) + a(t). Both trajectories must satisfy the 
equations of motion (3). If the variation a is assumed to 
be infinitcsimally small and the equations of motion are 
linearized in <j, we obtain the variational equations 
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FIG. 5. The ScPO (black) and its branches (green and brown) 
at E = -2.3 eV. 



where V is the Hessian matrix of the Hamiltonian, 

d 2 H 

%} dxidxj ' 

We integrate the combined systems (3) and (16) with 
arbitrary initial conditions Xq and ao, and we ask how 
fast the length of the tangent vector a(t) will grow. The 
Lyapunov exponent of the trajectory starting at xq is 
defined by 

A(. 0l ao)= ( limilnM i ( 17) 

where ||er|| denotes the length of the vector a. This defini- 
tion corresponds to an exponential growth ||c(t)|| oc e xt . 
In general, the tangent vector a will quickly align itself 
with the direction in which the expansion rate is largest. 
The resulting Lyapunov exponent does then not depend 
on the arbitrarily chosen initial vector a . There is an 
exception, however, for a trajectory in an invariant man- 
ifold: If the vector a is initially chosen tangent to the 
invariant manifold, it will remain tangent to it at all 
times. In this situation, we can meaningfully compute 
a Lyapunov exponent parallel to the invariant manifold 
and a Lyapunov exponent in the full phase space. The 
invariant manifold is normally hyperbolic if the latter is 
larger than the former. 

The Lyapunov exponents are particularly easy to com- 
pute for a periodic orbit 40 . Because the evolution equa- 
tion (16) is linear in the variation vector a, its solution 
can be written as a(t) = Y(t) ■ <Jo with a matrix Y(t) 
that does not depend on a. The matrix Y(T) is called 
the monodromy matrix of the corresponding periodic or- 
bit with period T, its eigenvalues m lr ..,m 2 „ are the 
Floquct multipliers. For a periodic orbit with period T, 
we have Y(jjT) = {Y{T)Y for fj, = 1, 2, . . . . Therefore 



& = -J-V-a, a (t ) — (T , (16) 



a(fiT) = (Y(T)r ■ <T . (18) 



FIG. 6. The BPO and its branches at E = 



-3.3 eV. 



So, m^, 



are the eigenvalues of Y(/zT). The spec- 



trum of Lyapunov exponents of the particular periodic 
orbit is then 



A,- = lim In Imf I = — 



In I m, 



(19) 



and the largest of the Floquet multipliers m, will give 
the Lyapunov exponent (17). For a periodic orbit in an 
invariant manifold we can use the eigenvectors of Y(T) to 
distinguish whether eigenvalues correspond to variations 
parallel or transverse to the invariant manifold, and we 
can then choose the largest Lyapunov exponents in the 
parallel and transverse directions. 

In a Hamiltonian dynamical system, the eigenvalues 
of the stability matrix Y(T) will always occur in pairs 
6 ±at Qr e ±i^T w jt n rea j numbers A and ip. These types 
of eigenvalues correspond to variations in unstable and 
marginally stable directions, and yield Lyapunov expo- 
nents A and 0, respectively. A third possibility arises in 
Hamiltonian systems with three or more degrees of free- 
dom: Eigenvalues can occur in quartets e (± A±i ¥>) T . This 
case is not relevant for our situation because the periodic 
orbits we are studying lie within the centre manifold, 
which is a subsystem with only two degrees of freedom. 

The parallel and perpendicular Lyapunov exponents 
will in general be different for different trajectories in the 
invariant manifold, though they will be equal for trajec- 
tories on the same invariant torus or in the same chaotic 
sea. To verify normal hyperbolicity numerically, we must 



TABLE I. Floquet multipliers of the Symmetric Stretch Peri- 
odic Orbit (SSPO). 



Energy 


Within centre manifold 


Off centre manifold 


-4.35 


-0.949855 ±0.312691i 


543.591 0.00183962 


-4.32547 


-1 ±0.000746992i 


513.284 0.00194824 


-4.02482 


-2.80253 


-0.35682 


2.80224 0.356857 


-4.02425 


-2.8049 


-0.356519 


0.950972 ± 0.309275i 


-4.02251 


-2.81218 


-0.355595 


-0.98557 ± 0.16925i 


-4.02195 


-2.81447 


-0.355306 


-2.83866 -0.352278 


-4.0 


-2.9049 


-0.344246 


-52.8575 -0.0189188 


-3.5 


-4.63902 


-0.215563 


-1922.55 -0.000520141 


-3.0 


-5.8626 


-0.170573 


-6579.25 -0.000151993 



therefore calculate Lyapunov exponents for a large num- 
ber of representative trajectories and check that the per- 
pendicular Lyapunov exponent is larger than the parallel 
exponent in all cases. 

Earlier studies of the dynamics in the collinear 
subsystem 7,10 found an energy interval in which the 
SSPO is stable against variations within that subsys- 
tem, which is transverse to the centre manifold. Also 
the narrow energy interval coincides with those values 
found recently by Inarrea et al 36 in the collinear case for 
the Porter-Karplus potential energy surface. For these 
energies, the transverse Lyapunov exponent of the SSPO 
is zero, and it is clear that this situation must violate the 
condition of normal hyperbolicity as soon as the trans- 
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(a) E = -4.35eV 




3.39 3.40 3.4] 3.42 



r (a.u.) 

(c) E=-4.3eV 




3.2 3.3 3.4 3.5 3.6 



r (a.u.) 



(e) E= -3.9 eV 




2.8 3.0 3.2 3.4 3.6 3.8 4.0 4.2 



r (a.u.) 



(b) E = -4.32547 eV 




3.0 3.2 3.4 3.6 3.8 4.0 



r (a.u.) 

(f) E=-2.3eV 

1.5 p" 




-1.5 |- 

2.5 3.0 3.5 4.0 4.5 5.0 5.5 

r (a.u.) 



FIG. 7. Poincare surface of section y — 0,p y > for dynamics within the centre manifold at the energies (a) — 4.35eV, (b) 
-4.32547 eV, (c) -4.3eV, (d) -4.0 eV, (e) -3.9eV (BPO is unstable) and (f) -2.3eV (ScPO is unstable). The main periodic 
orbits are labeled in (a), (c) and (e). 
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-0.3 -0.2 -0.1 0.0 0.1 0.2 0.3 



y(a.u.) 

FIG. 8. Poincare surface of section r = 3.40166 a.u. for dy- 
namics within the centre manifold at the energy — 4.3eV. 



verse Lyapunov exponent decreases below that within the 
centre manifold. Table I shows the Floquet multipliers 
of the SSPO within and transverse to the centre mani- 
fold. The SSPO is unstable within the centre manifold 
for those energies where it is stable in the collinear direc- 
tion. As the energy increases further, the SSPO is unsta- 
ble in both directions, but the instability in the collinear 
direction grows faster than that within the centre mani- 
fold, so that soon the SSPO does not violate the normal 
hyperbolicity of the central sphere any more. 

On its own, this observation does not allow us to con- 
clude that the centre manifold returns to being normally 
hyperbolic. It remains possible that normal hyperbolic- 
ity could be broken by any orbit other than the SSPO. To 
check this, we have calculated the Lyapunov exponents 
for a variety of orbits in the centre manifold over a range 
of energies up to the conical intersection ridge. It turns 
out that across the entire range of energies no orbit apart 
from the SSPO violates normal hyperbolicity. 

As an example of these calculations, Fig. 9(a) shows 
the Lyapunov exponents within and off the centre man- 
ifold for the energy -4.023 eV, at which the SSPO is 
stable in the transverse direction, and for orbits on the 
line p r = in the Poincare surface of section. This sec- 
tion includes the SSPO, ScPO, BPO and both regular 
and chaotic nonperiodic orbits. Because the SSPO forms 
the boundary of the surface of section, the two points 
with the highest and lowest admissible values of r cor- 
respond to the SSPO. The figure shows that normal hy- 
perbolicity fails for these points, but not for any other 
orbits. The enlargement in Fig. 9(b) confirms this con- 
clusion. Note that even for orbits arbitrarily close to the 
SSPO the transverse Lyapunov exponent is nonzero. Be- 
cause the SSPO is unstable under variations within the 
centre manifold, an orbit that starts arbitrarily close to 
the SSPO will quickly move away from it, and its long 
term behaviour will be entirely different from that of the 
SSPO. For this reason, the Lyapunov exponents can be 



discontinuous at the SSPO. 

Fig. 9(c) focuses on the Lyapunov exponents within 
the centre manifold. They are much smaller than the 
transverse Lyapunov exponents, and the difference be- 
tween trajectories on regular islands or in a chaotic sea 
can clearly be seen. For regular trajectories, we would 
expect these Lyapunov exponents to be zero. The nu- 
merical results show small, but finite values instead be- 
cause the Lyapunov exponents were obtained by solving 
the equations of motion for a finite time only, whereas the 
definition (17) requires the limit of infinitely long simula- 
tion time. If the actual simulation time is increased, the 
resulting Lyapunov exponents become even smaller. 

We have so far focused only on the question whether 
the central sphere is normally hyperbolic, i.e., whether 
the ratio k of the transverse to the parallel Lyapunov ex- 
ponents is larger than one. In fact, the precise value of 
this ratio is also relevant because the fundamental theo- 
rems about normally hyperbolic invariant manifolds 20,21 
guarantee that the stable and unstable manifolds of a 
NHIM exist and are differentiable at least k times. This 
result is important if normal form transformations are 
used to compute these manifolds, as they often have 
been 19 ' 23,26,29 ' 30 . The normal form will effectively rep- 
resent the invariant manifolds by Taylor series, which 
requires the existence of sufficiently high derivatives. Be- 
cause derivatives of order higher than k are not known to 
exist, the use of high order normal forms is questionable 
if the ratio k is low. 

Fig. 10 shows the ratio k of Lyapunov exponents for the 
SSPO, which is the orbit that potentially violates normal 
hyperbolicity, for energies from the saddle point up to 
the conical intersection ridge. The ratio is infinite just 
above the saddle point because the Lyapunov exponent 
within the centre manifold is zero. It decreases from there 
and reaches zero when the SSPO is stable. It then rises 
again and reaches a nearly constant value of k rs 5. As 
a consequence, we can expect the central sphere and its 
stable and unstable manifolds to be at least four times 
differentiable at all energies, except in a narrow range 
around the interval in which the SSPO is stable in the 
transverse direction. 



VI. CONCLUDING REMARKS 

We have shown that in the hydrogen exchange reac- 
tion the central sphere exists for all energies below the 
conical intersection ridge and that it will possess stable 
and unstable manifolds for all energies outside a small 
interval. The surprising fact that the normal hyperbolic- 
ity of the central sphere is restored after it has been lost 
implies that the phase space structure fundamental to 
TST, which consists of the central sphere and its associ- 
ated reaction tubes, will be in place even at energies high 
above the reaction threshold. Unfortunately, this result 
does not imply that the dynamics in the transition region 
will be simple. There will be homoclinic and heteroclinic 
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FIG. 9. (Color online) The Lyapunov exponents within centre manifold (blue) and off centre manifold (red) through the section 
p r = 0.0 in the Poincare surface for energy E — — 4.023 eV in the middle of the stable interval of SSPO. 



tangles that lead to complex phase space geometry and 
consequently to complex dynamics. In the collinear sub- 
system of the full system, this complex behaviour has 
been shown by Davis 13 and the most recently by Iharrea 
et al 36 . 

In the collinear subsystem, it is known that dynam- 
ics is as simple as assumed by TST only if the PODS is 
unique 7 . However, additional periodic orbits arise at en- 
ergies even lower than the energy at which the SSPO 
becomes stable, and trajectories that violate the no- 
recrossing assumption of TST appear at the same energy. 

Non-TST behavior in the full three-dimensional system 
must be at least as prevalent as in the two-dimensional 
subsystem. This means that even at energies at which 
the central sphere is normally hyperbolic, non-TST be- 
havior must be present. These energies are both below 
and above the range in which normal hyperbolicity is bro- 
ken. Thus, while the results of the current paper demon- 
strate that the normal hyperbolicity of the central sphere 
is more robust than one might have anticipated, this ro- 
bustness also implies that there is no direct link between 
the failure of TST and the violation of normal hyperbol- 
icity. It now becomes a separate question to determine 
what dynamical effects, and what phase space structures, 
cause the failure of TST. We will address this question 



in a forthcoming publication. 
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Appendix: The derivation of the kinetic energy expression 

The kinetic energy is derived from that given, for ex- 
ample, by Waalkens et al 26 . They study the HCN/CNH 
isomerization reaction in Jacobi coordinates: r the dis- 
tance between C and N, R the distance between H and 
the centre of mass of CN and 7 the angle between H, the 
centre of mass of CN and C (i.e.: The atoms C, N and 
H take the places of HI, H2 and H3, respectively, in our 



12 





FIG. 10. The ratio k of two Lyapunov exponents within and 
off the centre manifold for SSPO. 



Fig. 1). The corresponding kinetic energy expression is 



rr 1 2 



l 

2m 



i 

fir 2 



1 

mR 2 



(A.l) 



where /i = mcm^ / '(mc + mpf) is the reduced mass of 
CN and m = muivtic + THn)j(mH + mc + "ijv) is the 
reduced mass of the full system. In the exchange hydro- 
gen reaction, we have three identical atoms. Thus ft and 
m become \mn and |mjj, respectively. As a result the 
kinetic energy has the form 



rr 1 2 

T = p r 



4m 



H 



1 



p\ 



(A.2) 



It is singular when R = 0. This is the case for sym- 
metric collinear configurations such as the saddle point 
that is of central importance in our study. To avoid this 
singularity, we replace the polar coordinates R and 7 by 
Cartesian coordinates x and y, as shown in Fig. 1. The 
coordinate systems are related by 



R z 



i?COS7, 



V 



y = i?shi7, 
7 = arctan ^ 



We use r as the third coordinate as before. 

The generating function W associated with this trans- 
formation is 

W = p r r + p x R cos 7 + pyRsmj. 

It yields the following transformation of momenta: 
dW 



PR 



p 1 



dR 



dW 



p x cos7+p y sm7 

xp x + yp y 
R ' 



-p x Rsinj + pyRcosj 

= -ypx + xp y . 

Substituting these results into (A.2), we get 



T 



1 



Pr 



;(Px+Py) 



2n , ( x Py - VPxf 



which is the result used in (10). 
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